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Abstract. Low-luminosity AGNs (LLAGNs) are a very important class of sources since they occupy a significant 
fraction of local galaxies. Their spectra differ significantly from the canonical luminous AGNs, most notably by 
the absence of the "big blue bump". In the present paper, taking a typical LLAGN-NGC 4258-as an example, 
we investigate the origin of their spectral emission. The observational data of NGC 4258 is extremely abundant, 
including water maser emission, putting very strict constraints to its theoretical models. The infrared (IR) spec- 
trum is well described by a steep power-law form f v cx v~ , and may extend to the optical/UV band. Up until 
now there is no model which can explain such a steep spectrum, and we here propose a coupled jet plus accretion 
disk model for NGC 4258. The accretion disk is composed of an inner ADAF (or radiatively inefficient accretion 
flow) and an outer standard thin disk. A shock occurs when the accretion flow is ejected out of the ADAF to form 
the jet near the black hole, accelerating the electrons into a power-law energy distribution. The synchrotron and 
self-Comptonized emission from these electrons greatly dominates over the underlying accretion disk and can well 
explain the spectrum ranging from IR to X-ray bands. The further propagation of the shocked gas in the jet can 
explain the flat radio spectrum of NGC 4258. Several predictions of our model are presented for testing against 
future observations, and we briefly discuss the application of the model to other LLAGNs. 
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1. Introduction 

The nearby active galactic nucleus (AGN) in NGC 4258 
(M106) offers us a rare opportunity to understand the 
powering of a class of AGNs because of its precisely de- 
termined black hole mass and distance. High resolution 
22 GHz VLBI observations of the water maser emis- 
sion reveal a thin circumnuclear disk in Keplerian orbit 
(Greenhill et al. 1995a; Miyoshi et al. 1995; Moran et al. 
1995; Herrnstein, Greenhill, & Moran 1996). The maser 
emission extends from 0.13 pc to 0.26 pc, and the perfect 
Keplerian rotation curve requires a central binding mass, 
M, of (3.5 ± 0.1) x 10 7 M o within 0.13 pc. This makes it 
one of the two strongest supermassive black hole candi- 
dates to date (the other being our own Galactic center, 
see Melia & Falcke 2001 for a recent review). By measur- 
ing maser spot acceleration or proper motions (Greenhill 
et al. 1995b; Herrnstein et al. 1999), the distance of the 
source can also be directly determined independently of 
Hubble constant, and is found to be d = 7.2 ± 0.3 Mpc. 

In addition to the black hole mass and distance, the 
third critical parameter for AGN modeling is the mass ac- 
cretion rate M. The accretion rate in the maser disk can 
be determined from the observed water maser emission. 
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The projected positions of individual maser spots indi- 
cate that the maser disk is significantly warped (Miyoshi 
et al. 1995; Greenhill et al. 1995a). Such a warped disk 
will inevitably be illuminated obliquely by the central 
X-ray source, and Neufeld, Maloney, & Conger (1994) 
show that X-ray irradiation of dense molecular gas is 
an effective mechanism for generating powerful water 
maser emission. This also naturally explains the associ- 
ation of water megamasers with AGNs (Claussen & Lo 
1986; Braatz, Wilson, & Henkel 1994). By modeling the 
molecular annulus as a standard thin accretion disk (SSD; 
Shakura & Sunyaev 1973) illuminated by the X-ray pho- 
tons, Neufeld & Maloney (1995, hereafter NM95) infer 
the mass accretion rate through the disk to be Mnm95 = 
7x 1O _5 oM0 yr -1 , here a Ss 1 is the conventional viscosity 
parameter. 



However the NM95 accretion rate is obtained from 
modeling the masering disk, which is very far away from 
the hole since 0.13 pc ~ 40000r s with r s = 2r g = 
2GM/c 2 is the Schwarzschild radius of the black hole. 
The viscous timescale at 0.13 pc is extremely long, ~ 
10 9 al}(200K/T)(r/0.2 pc) 1 / 2 yr, therefore, it is very un- 
likely that accretion has been steady during such period 
(Gammie, Narayan, & Blandford 1999, hereafter GNB99). 
It is thus likely that the accretion rate at the inner region 
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of the accretion disk can has a different value from NM95 
(GNB99). 

There is a plethora of observational data for 
NGC 4258, ranging from the radio, IR, optical/UV, 
to X-ray bands, as we will review in Sect. 2. These 
broadband data indicate that NGC 4258 is very under- 
luminous, L/Le ~ 10~ 4 , where Le is the Eddington lu- 
minosity, making NGC 4258 an ideal candidate for an 
advection-dominated accretion flow (ADAF; see reviews 
by Narayan, Mahadevan, & Quataert 1999 and Kato, 
Fukuc & Mineshige 1998). In particular, since a thin disk 
is observed at the region out of ~ 0.13 pc, a truncated 
ADAF model is viable in which an ADAF connects at 
a transition radius r tr to an outer SSD. This would natu- 
rally explain why there is no maser emission inside 0.13 pc. 
Several transition mechanisms have been proposed (e.g., 
Meyer & Meyer-Hofmeister 1994; Honma 1996; Czerny, 
Rozanska, & Zycki 2000; Spruit & Deufcl 2002), but a 
reliable prediction for the value of r tr is still lacking. 

Lasota et al. (1996) assume that such a hybrid ac- 
cretion disk is responsible for the spectrum of NGC 4258. 
This model is later refined in GNB99. The X-ray spectrum 
of NGC 4258, which is hard to explain by a SSD alone, 
can be well explained by the ADAF. The IR spectrum of 
NGC4258 was not certain at the time when GNB99 was 
published, although GNB99 could well explain the IR flux 
which was available at that time. The IR spectrum in this 
model is produced by the outer SSD, therefore a very hard 
IR spectrum with the classical form of f v cx i/ 1 ' 3 is pre- 
dicted. The required accretion rate M ~ 0.01M Q yr _1 
(GNB99). This is more than two orders of magnitude 
higher than the NM95 rate. 

After the publication of GNB99, new IR observations 
determined a better spectrum (see next section for de- 
tails), which is in serious conflict with the prediction of 
GNB99. The detected spectrum is very steep, with spec- 
tral index s = 1.4 ± 0.1 (f v cx v~ s ). On the other hand, 
when fitting the spectrum of NGC4258 (or any source), an 
ADAF model usually does not take into account the role of 
the jet, which is believed to occur in the innermost region 
of the accretion disk. In some cases, the emission from the 
jet may dominate over the underlying accretion disk in 
sources such as blazars, M87 (Wilson & Yang 2002), and 
some classes of X-ray binaries (Markoff, Falcke, & Fender 
2001). 

In the present paper, we propose an alternative model 
to explain the spectrum of NGC 4258. We suggest that 
the main emission of NGC 4258 comes from the base of 
its jet rather than the accretion disk. In next section, 
we will present a review of the relevant observations of 
NGC 4258. Sect. 3 presents our model and explanations 
for the spectrum of NGC 4258, along with some theoreti- 
cal predictions. The application of our model to LLAGNs 
in general is discussed in the last section. 

2. Review of Observations 



2.1. Radio 

VLA observations reveal a compact radio source within 
0. 1 of the maser emission, with a ~ 3 mJy flux density 
at 2 cm, and a spectral index s of ~ —0.4 (Turner & Ho 
1994), indicated by the short "slanted" line in Fig. 1. 

Hcrrnstein et al. (1998) report that the 22 GHz radio 
emission traces a jet-like structure oriented along the ro- 
tation axis of the maser disk, and composed of distinct 
northern and southern components with flux density of 
2.5 — 3.5 mJy and 0.5 mJy, respectively. Because the rel- 
ative position of the center of mass of the molecular disk 
can be measured to a fraction of a milliarcsecond with 
VLBI, NGC 4258 provides a rare opportunity to circum- 
vent the ambiguity associated with core-jet emission. The 
centroids of the northern and southern components appear 
at 0.35 — 0.46 mas and —1.0 mas from the disk center, re- 
spectively. Both the flux and centroid of the northern com- 
ponent vary significantly with time. More interestingly, 
Hcrrnstein et al. (1998) report a 3 a upper limit of 220 /zJy 
on any 22 GHz continuum emission coincident with the 
central mass of the NGC 4258 maser disk. This result 
places an extremely strict constraint on theoretical mod- 
els of NGC 4258, and is the main motivation for GNB99's 
revisiting of the ADAF model after Lasota et al.'s (1996) 
work. GNB99 argue that such an upper limit requires the 
SSD- ADAF transition radius r tr ~ (10 - 100)r s . 

2.2. Optical/UV 

Wilkes et al. (1995, hereafter W95) have measured the 
5500 A flux toward NGC 4258 in polarized light. This 
emission presumably arises in the central engine and has 
been scattered into our line of sight. They find that the 
spectral index is well- fit by a power law with s = 1.1 ±0.2. 
This value agrees well with the argument of Stuwe, Schulz, 
& Hiihnermann (1992) that a continuum with s ~ 1.2 
can produce sufficient ionizing power to explain the ob- 
served strength of the optical emission lines in NGC 4258. 
The corresponding 5500 A central engine luminosity is 
highly sensitive to the type of scattering screen invoked. 
W95 proposed two candidates for the polarization: scat- 
tering by electrons in an H II region surrounding the nu- 
cleus, and scattering by dust grains located above and 
below the molecular disk. The estimated luminosities are 
£5500 ~ 10 39 ergs~ 1 A _1 and 10 37 ergs _1 A _1 , respectively. 
We denote both of them in Fig. 1 with empty circles be- 
cause they are not very precise. However, we think that 
the latter is more likely because connecting the near-IR 
data to the former point would result in an unbelievably 
hard spectrum, in serious conflict with the observed po- 
larization spectrum. 

2.3. Infrared 

Chary & Becklin (1997) present J, H, and K band Keck 
observations. Using the J band surface brightness profile 
as a baseline, they obtain the J-band subtracted H and K 
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bands image which is consistent with a unresolved point 
source. The nucleus is then found to be very red. Such 
redness can be attributed to either a red central source or 
dust extinction of starlight. They argue for the former by 
showing that this color cannot be explained by uniform 
foreground dust extinction, assuming a late-type stellar 
population typical of the bulge. Assuming that the intrin- 
sic spectrum has a similar power-law form with the op- 
tical/UV polarization spectrum obtained by W95, which 
gives f v oc ^- 11±0 - 2 ; they deduce the extinction and lu- 
minosity. 

Obviously such a spectral result, especially the value 
of the power-law spectral index, is not so certain. For ex- 
ample, a power-law fit of f v oc v 1 ^ also works well, as 
pointed out by GNB99. Using HST (NICMOS), the 5.1 
m Hale telescope and the Keck II telescope (MIRLIN), 
Chary et al. (2000) went on to present high-resolution im- 
ages from 1 to 18 /jm, placing a much better constraint 
on the spectrum. They find that the mid-IR data before 
extinction-correction share a similar spectrum with the 
optical/UV polarization results. Since mid-IR suffers rela- 
tively little dust extinction, such a similarity strongly sug- 
gests that the IR spectrum has approximately the same 
power-law form as the optical/UV polarization spectrum. 
To make sure, they tested four predictions for the 3.5 /im 
flux based on four different power-law spectral indexes, in- 
cluding s — —1/3. They find that only a steep power-law 
with a similar spectral index as the polarization spectrum 
of W95 can give a correct prediction with their obser- 
vation (R. Chary private communication). In addition, 
the extinction coefficient obtained from such a steep spec- 
tral fit across all the IR wavebands is consistent with the 
extinction derived from the N-band data point after in- 
cluding the 9.7 micron silicate feature (R. Chary, private 
communication). They therefore conclude that the intrin- 
sic extinction-corrected IR spectrum of NGC 4258 can be 
well fitted with a steep power-law form, and their best 
fit extinction-corrected spectrum goes as f v oc v~ i- 4 ^ - 1 . 
Considering W95, such a spectrum is very likely to extend 
to the optical/UV band. The data is listed in Table 1 and 
shown by filled circles in Fig. 1. 

Such IR and optical/UV spectral results, combining 
with the radio data, give an estimate the bolometric lumi- 
nosity of NGC 4258 of L bol « 8 x 10 41 - 3 x 10 42 erg s" 1 . 
Spinoglio & Malkan (1989) suggest that the 12 /zm lu- 
minosity is typically about 1/5 of the bolometric lumi- 
nosity of AGNs, independent of whether the emission 
is thermal or nonthermal in origin. This argument gives 
Lbol ~ 10 42 erg s^ 1 for the mid-IR data of Chary et al 
(2000). The consistency between the two estimates again 
suggests that the steep power-law spectrum extends to the 
optical/UV, and that there is no big-blue-bump as usually 
found in luminous AGNs. The steep power-law spectrum, 
/„ oc v~ 1A , found throughout the IR and possibly extend- 
ing to optical/UV strongly suggests a non-thermal origin, 
and that the dust in NGC 4258 is likely to be cool and 
contribute only at wavelengths longer than mid-IR (Chary 
et al. 2000). 



Such a steep IR spectrum is very unusual compared to 
the canonical luminous AGNs but seems to be common 
in LLAGNs, as we will illustrate in Sect. 4. It provides a 
strict constraint to the theoretical model of NGC 4258Q 
In particular, this is in serious conflict with the prediction 
of f v oc v 1 / 3 by GNB99, as discussed in Sect. 1. 

2.4. X-rays 

Several telescopes have observed the X-ray emission of 
NGC 4258 in recent years, from ASCA nearly eight years 
ago to XMM- Newton reported this year. The reader can 
refer to Table 1 in Pietsch & Read (2002) for the values 
of spectral parameters. An interesting result we find is 
that the observed X-ray spectral indexes differ greatly for 
different observations, ranging from s = 0.64 (Pietsch & 
Read 2002) to s = 1.11 (Fiore et al. 2001). We put these 
two most extreme results in Fig. 1. While one reason for 
the discrepancy of the spectral index is likely due to the 
instruments, it is possible that the spectrum is intrinsically 
variable with time. 

Reynolds, Nowak, & Maloney (2000) report their 
ASCA detection of the iron Ka line. They find the line 
is narrow^], which constrains the inner radius of the SSD 
to be > 100r g . Note in GNB99 the SSD-ADAF transition 
radius is r tr ~ (20 — 200)r 9 , so it only marginally satisfies 
the iron line observation. 

In the BeppoSAK observation of Fiore et al. (2001), 
large amplitude (100%) variability is observed on 
timescales of a few tens of thousands of seconds, while 
variability of ~ 20% is observed on timescales as short as 
1 hr. The size of the X-ray source is then constrained to 
be < 100r ff . 

2.5. Summary of Observations 

In short, the spectrum of NGC 4258 is quite different from 
the usual luminous AGNs. As discussed above, the data 
puts very strict constraints on the theoretical models and 
challenges our understanding of the powering mechanism 
in NGC 4258. On the other hand, we note that NGC 4258 
is somewhat analogous to Sgr A*: both sources have flat 
radio and steep IR spectra, and both are low luminos- 
ity sources with the luminosity peaked at sub-mm or IR 
rather than optical as in luminous AGNs. However, the 
bolometric luminosity of NGC 4258 is about 10 4 times 
higher than Sgr A* in units of Eddington luminosity, and 
the jet in NGC 4258 is much stronger than in Sgr A*, 

1 Chary et al. (2000) did try to fit such a spectrum using an 
ADAF model, but their calculation is too simplified. A more 
exact ADAF calculation predicts an exponentially decreasing 
spectrum at IR band, as shown by, e.g., the dot-dashed line in 
Fig. 1 in the present paper. 

2 The detected iron Ka line is moderately strong. This is 
somewhat unusual since among LLAGNs the iron line is either 
extremely weak or absent (Terashima et al. 2002). The iron 
line in NGC 4258 may come from, e.g., the warping molecular 
disk. 
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Table 1. Radio, infrared, and optical observational data for NGC 4258. 



Frequency (Hz) 


Flux(mJy) 


Reference 


— 10 X 1U 


o[S — -\J.Q ) 


i 

i 


2.2 x 10 10 


2.5-3.5, < 0.22'' 


2 


1.68 x 10 13 


435 


3 


2.40 x 10 13 


270 


3 


2.86 x 10 13 


195 


4 


8.70 x 10 13 


46 


3 


1.36 x 10 14 


25 


3 


1.88 x 10 14 


16 


3 


2.73 x 10 14 


10 


3 


5.45 x 10 14 


1.7,170 c 


5 



Notes and references: 

a : s=-0.4 is the spectral index (/„ oc v~ s ) 

b : 2.5-3.5mJy is the flux detected from the jet, while 0.22mJy is the upper limit detected from the location coincident with the 
central mass of the maser disk 

c : the two values correspond to two scattering models 

(1) Turner & Ho 1994; (2) Herrnstein et al. 1998; (3) Chary et al. 2000; (4) Rieke & Lebofsky 1978; (5) W95 



which if present has not yet been detected. The spectral 
similarity between the two sources does, however, sug- 
gest a common physical mechanism. A jet-ADAF model 
has been found to be very successful at explaining the 
spectrum of Sgr A* from radio to hard X-ray (see Yuan, 
Markoff, & Falcke 2002 for details; hereafter YMF02)[]. 
We consider that the same type of model may be applica- 
ble to the spectrum of NGC 4258 as well. In this model, 
most of the emission comes from the base of the jet rather 
than the accretion disk. The much higher luminosity of 
NGC 4258 compared to Sgr A* is then due to its more 
powerful jet. 



and the temperature of the flow will be increased nearly 
to the thermal temperature ~ 10 12 /(r/r g )K, since the en- 
ergy loss from the pre-shock ADAF is very small and can 
be neglected. Most of the emission from such a disk-jet 
system then stems from the high-temperature post-shock 
plasma. We find that the shock is radiative, i.e., the ra- 
diative cooling timescale is much shorter than the char- 
acteristic shock propagation time. We will show that the 
synchrotron emission will be responsible for the steep IR 
spectrum of NGC 4258 while its inverse Compton emission 
can explain the X-ray spectrum. The flat radio spectrum 
is explained by the emission from the outer part of the jet. 



3. Jet Model for NGC 4258 

The physical scenario of our model is as follows. The ac- 
cretion flow outside of 0.13 pc (or smaller) is described 
as a SSD. At a certain radius, a SSD-ADAF (or more 
generally a radiatively inefficient accretion flow) transi- 
tion occurs. Since the accretion onto the black hole must 
be transonic, the accretion flow becomes radially super- 
sonic after passing through a sonic point r Bon i c ^ 10r g 
(e.g., Narayan, Kato, & Honma 1997). When the radially 
supersonic accretion flow is subsequently ejected out of 
the disk to form the jet, a standing shock occurs due to 
the bending, as suggested by the general relativistic MHD 
numerical simulations of jet formation (Koide, Shibata & 
Kudoh 1998; Koide et al. 2000). After the shock, the ki- 
netic energy of the pre-shock accretion flow will be trans- 
ferred into the internal energy of the post-shock plasma 



Recently the presence of a four-month cycle in the radio 
variability of Sgr A* was reported (Zhao, Bower, & Goss 2001; 
see also Yuan & Zhao 2002). By analogy with the plateau state 
of GRS 1915+105, it is suggested that such quasi-periodic vari- 
ability is very likely to be produced through a jet (Zhao, Bower, 
& Goss 2001). If verified, this observation would strengthen ar- 
guments for a jet component in models of Sgr A* . 



3.1. The Shock at the Jet Base 

We denote the shock (or jet) radius as r$. Obviously, tq 
should be less than the sonic radius of ADAF, tq Si r son i c . 
Some accretion flow outside of r son i c may also go into the 
jet, but since it will not be shocked, its radiative contri- 
bution is small enough to be neglected. 

We calculate the temperature of the post-shock plasma 
from the parameters of the pre-shock flow using the shock 
transition conditions. The post-shock equivalent temper- 
ature is determined by 



rpeq 
1 2 



(r + i) m p 



1 m 
5~fc 



P -vl{r ). 



(1) 



2(r + l) 2 k 

The total thermal energy density of the post-shock plasma 
is then 



e = SnkT^ . 



(2) 



Here the strong shock approximation and a parallel shock 
configuration are assumed. The adiabatic index T — 1.5 
since the protons are nonrelativistic while the electrons are 
relativistic, n(= n p = n e ) is the particle number density, 
vi (ro) is the pre-shock velocity of the accretion flow at the 
shock location ro. We obtain Vi(ro) by solving the global 
solution for ADAFs. 
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However, since we do not know what fraction of the ac- 
cretion flow will go into the jet, and because hydrostatic 
equilibrium in the vertical direction of the disk does not 
hold at the jet formation region, it is hard to determine the 
density of the post-shock plasma. As in YMF02, here we 
adopt the mass loss rate in the jet, Mj ct , as a free param- 
eter. Then the density of the post-shock plasma follows 
from 



Mjot = iTTr n p m p V2- 



(3) 



Here V2 is the velocity of the post-shock flow and we calcu- 
late its value from the parameters of the pre-shock accre- 
tion flow using the shock transition conditions. Obviously 
Mjet should be smaller than the disk accretion rate M. 

We assume that electrons are accelerated in the shock 
to a power-law energy distribution, with a minimum 
Lorentz factor 7 m i n : 



iV(7e) = ^07e~ 



(je > 7min)- 



(4) 



Generally p « 2.5 is assumed in the literature for the in- 
jected distribution in the case of strong relativistic shock 
acceleration, e.g., Gamma-ray Bursts (GRBs). But in the 
present case, the shock is sub-relativistic and only moder- 
ately strong. For such a shock, we still know little about 
the form of the energy distribution of accelerated/heated 
electrons (e.g., Draine & McKee 1993). We conjecture that 
p is likely somewhat larger, or even a thermal distribu- 
tion is possible in the extreme case. We set p = 2.8 in 
the present paper to relate it to the observed IR spec- 
tral index s — 1.4 R As the result of radiative losses, the 
steady distribution becomes steeper, i.e., p — > p + 1 = 3.8 
(Kardashev, Kuz'min, & Syrovatskii 1962). 

The value of 7 m ; n is important for our spectral fitting. 
From Eq. (1), we can only determine the equivalent tem- 
perature of the post-shock plasma, but we do not know the 
respective temperature of electrons (or 7 m i n ) and protons. 
Another important quantity for determining the radiative 
cooling is the strength of the magnetic field. Both of these 
are difficult to estimate from first principles. We therefore 
follow the usual approach, using two parameters, e e and 
6b , to describe respectively the magnitude of the energy 
densities of random relativistic electron motions and of the 
magnetic field in terms of their ratios to the total thermal 
energy density of the post-shock plasma, 



P - 1 7min"- e ra e <^ 
p-2 e 



D 1 



(5) 



(6) 



4 In the case of Sgr A*, the IR upper limits require that p 
must be much larger if assuming a power-law form, or, the 
distribution is approximately a relativistic Maxwellian form 
(YMF02). The discrepancy between NGC 4258 and Sgr A* 
may be due to the different physical parameters at the shock 
region, especially the shock strength. On the other hand, due to 
the poor constraints on the IR spectrum in Sgr A*, we cannot 
exclude that some fraction of electrons may be accelerated into 
a very steep power-law distribution. 



To calculate the synchrotron and self-Comptonized 
emission from the post-shock plasma, we still need to know 
the profile of the energy density after the shock front. We 
find to fit the spectrum, the shock must be radiative. This 
means e e ~ 1 and so the radiative cooling time of electrons 
due to synchrotron emission, 



3 8tt m e c 

^rad — T 



(7) 



AcTTje^B 2 V100/ V 1000G / 
is much shorter than the dynamical time of the radiative 



region, 



tdy 



ro 

V2 



10' 



GST 



(8) 



for typical values of 7 e ~ 7min, B and V2 in our model. We 
find that this holds true even when the self-absorption of 
synchrotron emission is taken into account. 

In contrast to NGC 4258, the shock in Sgr A* is adia- 
batic. Besides the difference in mass between the two black 
holes, which will affect the dynamical timescale, this is 
due to the larger magnetic field in NGC 4258 compared 
to Sgr A*. This seems consistent with the idea of jets hav- 
ing a magnetic origin, because a stronger magnetic field in 
NGC 4258 is consistent with its jet being much stronger. 

The structure of radiative shocks has been discussed 
in many papers (e.g., Blandford & McKee 1976; Draine 
& McKee 1993; Granot & Konigl 2001). In addition to 
a "shock transition zone" where the bulk of the kinetic 
energy of the pre-shock gas is dissipated and which can 
be treated as an infinitesimal front, there exists a thin 
"radiative zone" or "cooling layer" after the shock front, 
where most of the dissipated energy of electrons is radiated 
away. In particular, the magnetic field can be amplified 
due to the cooling compression of the gas. Thus if the pre- 
shock magnetic field in the plasma is already sufficiently 
strong, as in our present case of ADAF, the magnetic field 
in the cooling layer can reach or even exceed equipartition 
with the electron energy density, i.e., e.e e e (e.g., Granot 
& Konigl 2001)0. 

Such a cooling layer is typically removed from the 
shock front by a short distance. An exact approach would 
be to calculate the radiation hydrodynamics of the post- 
shock plasma to obtain the distance-dependent profiles of 
n e , T e , and B after the shock front. We here simply treat 
the cooling layer as a very thin homogeneous slab imme- 
diately after the shock front. The emission from the jet 
downstream of this slab can be neglected since almost all 
the energy of electrons will be radiated away in this thin 
layer until the electrons are accelerated again (see Sect. 
3.3 below). The width of the cooling layer is determined 
by the post-shock velocity and cooling time of electrons, 



l c w u 2 irad = v 2 Sl c e e e/Lt t < r , 



(9) 



5 Such an approximate equipartition is also required in the 
modeling of blazars (Readhead 1994; Bower & Backer 1998), 
microquasars (Kaiser, Sunyaev, & Spruit 2000) and GRBs 
(Granot, Piran, & Sari 1999). 



6 



Feng Yuan et al.: NGC 4258: a jet-dominated low-luminosity AGN? 



where S is the "cross section" of the cooling layer. The 
effect of the vertical component of velocity of post-shock 
plasma is neglected since it should be small at the begin- 
ning of the jet. The shock may be oblique and the shock 
normal direction may be between the radial and vertical 
directions, and the cross section of the shock front is likely 
a ring rather than a full circle (thus the jet may have a co- 
coon structure). We here approximate the cross section of 
the shock as S — nr^. L to t is the emitted total luminosity 
from electrons, and is obviously a function of l c , which we 
therefore calculate iteratively. Note that we cannot simply 
use the ratio between the energy and synchrotron power 
to calculate the electrons' radiative lifetime i ra d because 
the effect of synchrotron self-absorption is significant. It is 
not, however, enough of an effect to change the "radiative" 
feature of the shock. 

For a given set of parameters n, e e and es, we numer- 
ically calculate the emitted luminosity and spectrum in- 
stead of using the simplified analytical formula because 
of the cutoffs in the emitting particle distributions. This 
is because the minimum frequencies of synchrotron and 

inverse Compton emission, t / syn(0 / niin) — 7min^ e ^/'^^"^^ e ^' 

and 7^ lin i / syn(7min), fall in the mid-IR and soft X-ray 
bands, respectively, and the slope of the inverse Compton 
emission will thus be changing in the fitted region, as dis- 
cussed in the next section. 

We would like to note that the number of free param- 
eters in our model is almost the same as the SSD-ADAF 
model of GNB99. There are five free parameters in our 
model, namely the shock radius r , the mass loss rate in 
the jet Mj et , the electron energy spectral index p, and the 
electron/magnetic energy parameters e e /es- As a com- 
parison, five free parameters are also required in GNB99, 
namely the accretion rate M, the SSD-ADAF transition 
radius tv, the viscous parameter a, the magnetic parame- 
ter /?, and the parameter describing the fraction of viscous 
energy that goes into the electrons, S. 

3.2. Infrared to X-ray Spectral Fitting Results 

As we stated in Sect. 2, the X-ray spectral slope varies 
between the different observations, e.g., s — 0.64 ± 0.08 
(Pietsch & Read 2002), 0.78 ± 0.29 (Makishima et al. 
1994), 0.87 ±0.15 (Reynolds et al. 2000), and 1.11 ±0.14 
(Fiore et al. 2001). In addition to instrumental reasons, 
it is very possible that such variations are intrinsic, as 
we know happens in some sources (e.g., Giebels et al. 
2002). We will therefore fit both the steepest (s = 1.11, 
from BeppoSAX) and the flattest (s = 0.64, from XMM- 
Newton) X-ray data in the framework of our model. 

Fig. 1 shows our best fit to the BeppoSAX results. 
The parameters are Mj ct — 2.9 x 1O~ 4 M yr _1 , r = 5r g , 
e e = 0.27 and es = 0.31. From these parameters we ob- 
tain B — 1400G, 7 m in = 30, the length of the cooling layer 
l c = 0.01r 9 , and the particle density n = 1.8 x 10 9 cm~ 3 . 
The long-dashed line in the figure shows the synchrotron 
emission, the dotted line shows its sclf-Comptonizcd radi- 
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Fig. 1. Model fit to the spectrum of NGC 4258. 
Parameters are chosen to fit the BeppoSAX data, which 
is described in Sect. 2. The thick solid line shows the 
emission from the radiative shock at the base of the jet, 
with the long-dashed and dotted lines being for the syn- 
chrotron and its self-Comptonized emission, respectively 
(Sect. 3.1 & 3.2). The thin solid line shows the emission 
from the outer part of the jet (Sect. 3.3). The dot-dashed 
and short-dashed lines are for the underlying ADAF and 
SSD, respectively, with two different accretion rates, and 
are minor contributors to the spectrum. 

ation, and the thick solid line shows their sum. The power- 
law IR spectrum with s = 1.4 is fitted by the synchrotron 
emission from electrons with a steady energy distribution 
of index p = 3.8 = 2s + 1. The predicted flux at the opti- 
cal band is also roughly in agreement with the lower flux 
estimated in W95, which we feel is more likely (see Sect. 
2.3). In the radio band, due to the strong synchrotron self- 
absorption and the fact that the emission from the region 
out of the cooling layer can be neglected, the synchrotron 
spectrum steeply decreases below the synchrotron peak so 
the 220 [iJy upper limit at 22 GHz is easily satisfied. 

The X-ray spectrum is mainly due to the synchrotron 
sclf-Compton component. Due to the effect of the low- 
energy cutoff in the electron distribution, the power-law 
spectrum at 2 — 10 keV is harder than s = (p— 1)/2 = 1.4. 
The exact value of the spectral index depends on the min- 
imum frequency of inverse Compton emission, ~ ^lynTmin' 
here ^f yn is the peak frequency of the synchrotron emis- 
sion, which is mainly determined by the value of B at the 
cooling layer. If ^ yn 7mi n i s smaller, the spectrum would 
become steeper. 

Also shown in Fig. 1 is the emitted spectrum from the 
underlying accretion disk, namely the ADAF (dot-dashed 
lines) and SSD (short-dashed lines). Because the value of 
the accretion rate is uncertain, we here adopt two values: 
M = Mnm95 (for the lower set of lines), and 10Mnm95 
(for the upper set of lines). The transition radius is set 
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Fig. 2. Model fit to the spectrum of NGC 4258. The pa- 
rameters are chosen to fit the XMM-jVeioion result. The 
lines are defined as in Fig. 1. The discrepancy in the mid- 
IR may be due to variability (see text for details). 



to be rt r =0.13pc, but our results are not sensitive to its 
exact value. The spectrum of ADAF is calculated by self- 
consistently solving the radiation hydrodynamic equations 
(e.g. Yuan et al. 2000). The ADAF with M = M N M95 
produces very low radio flux at 22 GHz because of small 
accretion rate, well satisfying the 220 /iJy upper limit. The 
accretion rate 10Mnm95 is the largest rate for which an 
ADAF can satisfy the 22 GHz upper limit for r tr = 0.13 pc 
0. In both cases, the contributions of ADAF and SSD to 
the spectrum of interest are rather small and can even be 
neglected. Decreasing rt r will increase the contribution of 
SSD, but we find even for rtr ~ 50r g — which is prohibited 
by the iron line observation — the contribution of SSD is 
still too small to contribute to the spectrum for M ^ 

10M N M95- 

Fig. 2 shows our results of fitting the XMM-Newton 
X-ray data. The parameters are r — 3r s , M- Jct = 7 x 
1O _5 M yr _1 , e e = 0.2, and e# = 0.11. Such parameters 
give B = 1250G, 7 min = 75, n = 1.2 x 10 9 cm" 3 , and 
l c = 0.004r g . The harder inverse Compton spectrum com- 
pared to Fig. 1 is due to the larger Vgynlmin- Thus, we in- 
terpret the variations in the X-ray spectrum of NGC 4258 
as due mainly to variations of 7min, which could result 
physically from changes in the shock location ro. However, 
compared to Fig. 1, this set of parameters significantly un- 
derpredicts the mid-IR flux of Chary et al. (2000). This 
may be due to intrinsic variability in NGC 4258, since 
the broadband data is not simultaneous. According to 
our model, when the X-ray spectrum becomes harder, 
the mid-IR flux is decreased due to the increased self- 
absorption from the required lower-energy electrons. Our 



results would therefore be consistent if the high mid-IR 
flux detected by Chary et al. (2000) was obtained when 
the X-ray spectrum was as steep as the BeppoSAK result. 
Future simultaneous IR and X-ray observations would be 
a good test of this model. 

The X-ray emission in our model comes from small 
spatial scales, ro ~ 3 — 5r g , and therefore rapid variability 
as detected in Fiore et al . (2001) would be expected. On 
the other hand, we can easily set r tr > 100r g so our model 
will not produce a broad iron line, consistent with the 
observation of Reynolds, Nowak, & Maloney (2000). 

In our model presented in Fig. 1., Mj et « 4Mnm95- If 
the accretion rate in the inner ADAF really equals Mnm95> 
this result is unphysical. It is possible that this factor of 
~ 4 discrepancy could be absorbed by the uncertainties 
in both the NM95 model and our own[|. However, as we 
show above, we can still be under the 22 GHz constraint 
with an accretion rate as high as 10Mnm95- In this case, 
a fraction of 25% (for Fig. 1) and 10% (for Fig. 2) of 
accretion flow will be transferred into the jet. We speculate 
that the latter is more likely. 

We now discuss why we require the accretion flow in- 
side r tr to be an ADAF (or more generally a radiatively 
inefficient accretion flow) rather than a SSD. To explain 
the spectrum of NGC 4258, especially the power-law IR 
spectrum whose luminosity is almost equal to the bolomet- 
ric luminosity of NGC 4258, extremely relativistic elec- 
trons with power-law energy distribution are required. An 
ADAF is basically adiabatic, i.e., the potential energy of 
the accretion flow is stored in the protons without being 
radiated away. When the flow is shocked at ro, the stored 
energy in the pre-shock ADAF will be transferred into 
the thermal energy of post-shock plasma. Among them a 
large fraction, e e w 0.2 — 0.3, will be used to accelerate 
electrons into a relativistic power-law distribution. Note 
that the radiative efficiency of the radiative shock is much 
higher than an ADAF, which follows, 



V 



Lboi 
M ict c 2 



,GMM\ 



jet 



(r - 2r g )M ict c 2 r /r g 



0.1-0.2.(10) 



The accretion rate can be even higher if the transition ra- 
dius is smaller, as in GNB99. 



Here we simply use the Paczyhski & Wiita (1980) poten- 
tial to mimic the geometry of the central black hole. If the 
ADAF were replaced by a SSD, the accretion energy would 
be radiated through a large range of radii. Therefore, even 
though a shock would still occur, only a small fraction of 
accretion energy would go into the thermal energy of the 
post-shock plasma. The electrons would not obtain a rela- 
tivistic power-law energy distribution and their radiative 
efficiency would be <C 0.1. 

7 For example, in addition to other simplifying approxima- 
tions in our model, we assume that the proton density is the 
same with the electron density, n p — n e when calculating n e 
from Mj ct . It is also possible that proton-proton collisions in 
the post-shock plasma, and the ensuing pion creation and de- 
cays, would lead to significant amounts of additional e . In this 
case n p < n e , and the required Mjet would become smaller. 
We are investigating this possibility elsewhere (Markoff et al. 
in preparation). 
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There exists another possible piece of evidence for the 
existence of an ADAF rather than a SSD within 0.13 pc. 
Observation indicated that within this radius the water 
maser emission disappears. In order to account for the 
absence of masers, NM95 propose that X-ray irradiation 
"shuts off" inside this radius. If this is true, Menou & 
Quataert (2001) show that for an accretion rate as low as 
A/nm95 , both MHD turbulence and gravitational instabil- 
ity will cease to be a viable angular momentum trans- 
port mechanism, therefore, the accretion in NGC 4258 
should proceed very inefficiently or perhaps not at all. 
This excludes the possibilities of producing the bolomet- 
ric luminosity of NGC 4258 through a SSD. The absence 
of masers may be due to the transition between SSD and 
ADAF. Thus the X-rays will irradiate both the SSD and 
the ADAF, and MHD turbulence will still be an efficient 
angular momentum transport mechanism even though the 
accretion rate in NGC 4258 is so low (Menou & Quataert 
2001). 



3.3. Radiation from the Outer Part of the Jet and 
Radio Observations 



According to the above scenario, a large fraction of accre- 
tion flow will be ejected out of the ADAF and go into the 
jet after passing through a radiative shock. The emission 
from the thin cooling layer downstream of the shock front 
explains the IR and X-ray spectra of NGC 4258. We now 
follow the propagation of the plasma and discuss the emis- 
sion from the jet downstream of the cooling layer. We con- 
sider only the northern jet in NGC 4258 since the southern 
jet is severely affected by the maser disk. If our model is 
correct, the propagated jet should be able to fit the ra- 
dio spectrum detected by Turner & Ho (1994), Herrnstcin 
et al. (1998), and the 22 GHz offset. We explain this by 
the canonical jet model of Blandford & Konigl (1979) (see 
also Falcke & Biermann 1999). The numerical approach is 
presented in Falcke (1996) and Falcke & Markoff (2000). 
The plasma will become supersonic after being acceler- 
ated in a "nozzle". The velocity profile along the jet is 
then calculated by solving the relativistic Euler equation. 
The opening angle of the jet is determined by the Mach 
number of the jet, ip ~ M~ l . Given the mass loss rate in 
the jet determined at the base of the jet, Mj et , the density 
profile is then determined. After the radiative zone, the 
jet plasma initially emits very little since the electrons 
are very cool after radiating most of their energy away. 
However, we assume that after propagating a certain dis- 
tance from the nucleus, zq, the electrons in the jet will be 
accelerated through, e.g., the internal shocks due to the 
velocity irregularities in the jet beam (Rees 1978). The 
electrons will then be accelerated again into a power-law 
form, N e j(je tj ) = N° tj j~?, with j e< j > 7 jjmin and tak- 
ing the standard strong shock value p = 2.5, although the 
model is not sensitive to the exact value. The average en- 



ergy of electrons (or equivalently 7j, m in) can be calculated 
by the relativistic shock transition condition, 

p — 1 

—r;l],mmm e c 2 = e j!e (Tj(z ) - l)m p c 2 . (11) 
P * 

Here Tj(z ) is the bulk Lorentz factor of the jet at Zq and 
6j e is the energy density factor of electrons. 

The shock location zq must be 3> r g . Herrnstein et 
al. (1998) obtained the 22 GHz upper limit by extracting 
the emission located within an ~ 0.24 x 0.24 mas box 
centered on the peak of the northern jet emission located 
at ~ 0.4 mas. This requires that the contribution to the 
22 GHz flux from the jet within the region < z < (0.4 — 
0.24/2) mas ~ 5600r ff must be below the 220/zJy upper 
limit. We therefore set zq — 6000r 9 , but again the fit is 
not sensitive to its precise value. 

Our calculation indicates that Fj(6000r g ) w 3 and 
7j,min ~ 650 for tj e = 0.5. Further assuming equipar- 
tition between the electrons and magnetic energy, and a 
viewing angle of 9 = 83° (Miyoshi et al. 1995), we can cal- 
culate the emergent spectrum of the jet and the location 
of the peak, z(v), of the emission at a frequency v. The 
result is shown by the thin solid line in Fig. 1. We see that 
it fits the results of Turner & Ho (1994) and Herrnstein 
et al. (1998) very well [| Moreover, in our result, the 22 
GHz emission mainly comes from the location ~ 0.4 mas 
from the nucleus, satisfying the observed (time variable) 
0.35 — 0.46 mas offset (Herrnstein et al. 1998). We explain 
the variation of the offset with time as due to the variation 
of the shock location zq. If zq <C 6000r s , to fit both the 
22GHz flux and its offset simultaneously, a much smaller 
Afjct would be required. But in this case, the synchrotron 
emission from the jet within the region < z < zq would 
greatly exceed the 0.22 mJy upper limit. This supports 
our adopted Mj ct , and further, our jet-ADAF model of 
NGC 4258. 

As shown by the thin solid line in Fig. 1, the contri- 
bution from the outer part of the jet to the X-ray flux 
can be neglected compared to the cooling layer at the 
base of the jet. NGC 4258 is the only source for which, 
to our knowledge, both the offset and an upper limit of 
the 22 GHz emission from the nucleus are so precisely 
determined. This suggests that it may be possible that 
the X-ray spectra of some sources which were previously 
attributed to emission from the outer part of the jet ac- 
tually comes from the base of the jet (or in other words, 
the coupling region between the jet and accretion disk). 
Of course, an additional important factor when we com- 
pare the respective contribution of the outer jet and base 
of jet is the Doppler factor S (5 = [Tj(l — f3jCos9)}~ 1 ). At 
the base of the jet, S « 1, but for the outer part of the 

8 We find the synchrotron emission in the jet in NGC 4258 is 
optically thin. This can be understood directly from observa- 
tions (GNB99). Observation indicated that the opening angle 
of the jet is about 10° - 15° (Herrnstein et al. 1997). The ra- 
dio brightness temperature is ~ 10 9 K, which is much smaller 
than the kinetic temperature of electrons responsible for the 
emission, T e ~ 7j,min»7ieC 2 /3A; w 10 12 K. 
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jet, S in general deviates from 1 significantly, depending 
on the value of 6. For NGC 4258, 9 » 83°, so S < 1, and 
the emission from the base of the jet dominates over that 
from the outer part of the jet. But for some sources like 
Blazars, (5 > 1, thus the emission from the outer part of 
the jet may dominate over that from the base of the jet. 
Taking the effect of <5 into account, the candidates in which 
emission from the base of the jet may dominate over the 
outer jet are those having misaligned jets. M87 will be a 
good example. 

3.4. Model predictions 

Taking into account the emission from the outer part of 
the jet and the base of the jet, our model makes the fol- 
lowing predictions for NGC 4258: 

— There exists an infrared bump peaked at ~ 50/im. The 
radio spectrum should be flat with spectral index of 
s k 0.3 until a break frequency k, 200 — 400GHz, 
with the exact value depending on the accretion rate 
of the underlying ADAF. These results can be tested 
by the future ALMA telescope. 

— Rapid variability at IR band with timescale of ~ 1 hr 
is expected as detected in X-ray band, and the IR and 
X-ray variability should be tightly correlated. 

— A somewhat weaker prediction of our model is that, 
if the mass ejection from the accretion disk suffers 
the same kind of ejection-related instability as in 
microquasar GRS 1915+105 ( Dhawan, Mirabel, & 
Rodriguez 2000) and maybe also Sgr A* (Zhao et al. 
2001), a similar quasi-periodic variability at submm- 
IR bands should be detected. The period in the case 
of Sgr A* is about four months (Zhao et al. 2001; Yuan 
& Zhao 2002). Scaled with the mass of the black holes, 
we predict a variability period of about 4^5 yrs at 
submm-IR bands in NGC 4258 . 

— Extending these predictions to other sources, we will 
argue in the next section that the spectra of some 
LLAGNs may have the same origin as NGC 4258. In 
this case we would predict strong correlations between 
the radio, submm, IR, and X-ray emission in these 
sources. 

4. Discussion: is NGC 4258 a typical LLAGN? 

Recently much attention has been paid to LLAGNs, a 
class which includes many of our nearby galactic neigh- 
bors. This class of AGN is very common in the local 
universe, e.g., 43% of all northern galaxies brighter than 
Bt — 12.5 mag are active in the form of emission-line 
nuclei classified as Seyferts or LINERs (Ho, Filippenko & 
Sargent 1997; Ho et al. 2001). The spectral energy distri- 
bution of LLAGNs are markedly different with the canon- 
ical luminous AGNs which are usually associated with the 
big blue bump (e.g., Elvis et al. 1994). Ho (1999) presented 
the spectra of a sample of seven LLAGNs for which rea- 
sonably secure black hole masses have been determined by 



dynamical measurements. Several most obvious features of 
their spectra are as follows (Ho 1999). 

— The optical/UV slope is quite steep. The averaged 
power-law index is 1.5, if two possibly highly reddened 
objects are excluded, while in canonical AGNs it is 
0.5-1.0. 

— The UV band is exceptionally dim relative to the opti- 
cal and X-ray bands. There is no evidence for a big blue 
bump component in any of the objects. The spectrum 
reaches a local minimum somewhere in the far-UV or 
extreme-UV region. 

— If fitting the X-ray spectra with a power-law (/„ oc 
v~ s \ the spectral index s obtained is very variable 
among the sources, s ~ 0.6 — 1.2 Q 

— There is tentative evidence for a maximum in the spec- 
trum at mid-IR or longer wavelengths. 

— The nuclei have radio spectra that are either flat or 
inverted. 

Ho (1999) argued that the absence of the big blue 
bump is a property intrinsic to the nuclei and not an ar- 
tifact of strong dust extinction. In addition to the above 
spectral features, the iron line in LLAGNs is in general 
extremely weak or absent (Terashima et al. 2002). 

In contrast with the obvious importance of LLAGNs, 
we are still not clear about their physics. From their low 
luminosity in units of Eddington luminosity, the absence 
of big blue bump, the weak or absent iron line, and other 
observational facts, we can deduce that the inner region 
of the accretion disk may be best described as an ADAF 
rather than SSD (e.g., Ho 2002). However, it is very hard 
for an ADAF to fit the spectra of LLAGNs. Quataert et al. 
(1999) use an SSD- ADAF transition model to fit the spec- 
tra of two LLAGNs, namely M81 and NGC4579. In ad- 
dition to underpredicting the radio spectrum, which they 
attribute to the "contamination by jets", the predicted 
optical/UV spectra are somewhat too steep. This is be- 
cause in their model the optical/UV is mainly produced 
by a SSD whose spectrum is exponentially decreasing at 
this waveband. It may be even more challenging to fit the 
spectra of other LLAGNs in the sample of Ho (1999) with 
such a model. 

On the other hand, the excellent agreement between 
the spectra of NGC 4258 and of general LLAGNs indi- 
cates that NGC 4258 is likely a typical LLAGN. In fact, it 
seems to be the LLAGN for which we have the best data 
so far, although Sgr A* is a close second. NGC 4258 is 
better, however, since in addition to the precise determi- 
nation of the mass of the black hole, the distance to the 
source, and the broadband spectrum ranging from radio 
to hard X-ray, NGC 4258 has additional extremely im- 
portant constraints. Specifically, the 220 /iJy upper limit 

9 After the publication of Ho (1999) paper, new X-ray obser- 
vational data of two sources are available, namely s = 1.2 for 
M87 (Bohringer et al. 2001; Wilson & Yang 2002), s = 0.88 for 
NGC4579 (Eracleous et al. 2002), and s = 0.94 for M81 (Page 
et al. 2002). 
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and the offset at 22GHz may provide ways to determine 
the mass accretion rate from the modeling of the water 
maser emission. Thus the success of our jet-ADAF model 
in explaining the spectra of both NGC 4258 and Sgr A* 
may be relevant to our understanding of other LLAGNs. 

We therefore suggest that the spectra of LLAGNs are 
at least partly due to the emission of the jet, namely that 
the synchrotron emission from the base of the jet is re- 
sponsible for the steep IR/optical/UV spectrum peaked 
roughly at mid-IR, and its self-Comptonized emission pro- 
duces the X-ray spectrum. The flat/inverted radio spec- 
trum comes from further out in the jet. In the radio wave- 
bands, this seems to be already well-established (Falcke 
2001; Ulvestad & Ho 2001; Nagar et al 2001). Several im- 
portant results of radio observations, i.e., high brightness 
temperature, the flat spectra, and the compactness of the 
radio cores, all are best explained by the emission from 
the jet. In the X-ray wavebands, the variety of the X- 
ray power law spectral indexes can be explained as due 
to the variety of both the peak frequency of synchrotron 
emission-^Pyjj-and Tmin among the sources, similar to the 
variations in the X-ray spectrum of NGC 4258 as we dis- 
cussed in Sect. 3.2. The variety of optical/UV power-law 
indexes may be the reflection of the complexity of shock 
type (radiative or adiabatic) and shock acceleration (the 
value of p) . 

The uncertainty in the extent of the underlying disk 
contribution adds complexity to this picture. It is possible 
that, at some wavebands or for some sources the emission 
from the underlying ADAF, and maybe also SSD, will be 
comparable with or even dominate over that from the jet. 
For example, observations indicate that X-ray variability 
on timescales less than a day is lacking in many LLAGNs 
although short-term variability has been detected in some 
LLAGNs (Ptak et al. 1998; Terashima et al. 2002). This 
may be an evidence for the dominance of ADAFs over 
jets at X-ray band in these sources because in an ADAF 
the X-ray emission can be dominated by bremsstrahlung 
process whose corresponding variability timescale is very 
long. The relative contribution of ADAF and SSD depends 
on the transition radius between them, while the relative 
contribution of the jet and accretion disk depends on our 
viewing angle with respect to the jet axis and the fraction 
of the accretion flow transferred into the jet. 
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